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In many industrial processes, such as pouring a liquid 
or coating a rotating cylinder, air bubbles are entrapped in- 
side the liquid. We propose a novel mechanism for this phe- 
nomenon, based on the instability of cusp singularities that 
generically form on free surfaces. The air being drawn into 
the narrow space inside the cusp destroys its stationary shape 
when the walls of the cusp come too close. Instead, a sheet 
emanates from the cusp's tip, through which air is entrained. 
Our analytical theory of this instability is confirmed by ex- 
perimental observation and quantitative comparison with nu- 
merical simulations of the flow equations. 




FIG. 1. a Cross section of the stationary air-fluid interface 
produced by a thin (1mm) stream of viscous oil poured into 
a deep pool of the same fluid. The position of the cusp is 
marked by a circle, b A hollow cylinder of air forms after 
the cusp has become unstable at a slightly higher flow rate. 
(Photograph by Itai Cohen) 

Air bubbles are a ubiquitous presence in fluid flow, ap- 



pearing when pouring a liquid into a beaker, when beat- 
ing an egg, or in river streams. This aeration is often 
desirable, for example to promote chemical reactions [0, 
yet in many industrial processes entrainment of air bub- 
bles is detrimental to the quality of the product, render- 
ing the flow unsteady. For example, it is the single most 
important factor limiting the speed at which paints or 
coatings can be applied to a solid surface dJ^. But in 
spite of its importance, no general understanding of air 
entrainment exists, except for the rather special circum- 
stance that the free surface conspires to enclose an air 
bubble from all sides, as was recently found for a dis- 
turbed water jet [Q. Instead, here we argue that the 
presence of cusp singularities on the free surface results 
in a generic mechanism for air entrainment. 

In recent years it has emerged that free surfaces are 
rather susceptible to the formation of sharp tips or cusps 
l^'iH^]- This is true in particular for viscous fluids, where 
shear stresses are strong compared with surface tension 
forces j7j. Examples are drop impact on a surface [ p^ , 
jets impinging on a pool of liquid Q], and the coating 
of a pre- wetted solid cylinder j^. Once the air enters 
these narrow passages, the gas flow serves to destroy the 
original structure to let the air penetrate below the sur- 
face. What is surprising about this novel mechanism is 
that the forces the gas flow exerts on the fluid plays the 
crucial role, even though the viscosity of the air may 
be smaller by many orders of magnitude than that of 
the fluid. Other instabilities may occur at a three-phase 
boundary, for example when the solid to be coated is dry, 
a problem studied in [^Ijsj. In the presence of surfactant, 
still another mechanism for the instability of a cusp has 
been suggested by Siegel 

For purposes of illustration, consider the particular ex- 
ample of a two-dimensional cusp that forms when a thin 
stream of a viscous silicone oil is poured into a container 
of the same fluid. Since the falling liquid drags other 
fluid away from the surface, a dip is produced around 
the fluid stream. Increasing the flow rate above a critical 
value, this dip is no longer smooth, but a singular point 
on the surface is approached with two vertical tangents. 
A cross section of this cusped profile is shown as a black 
silhouette in Fig.|l|a, the outer wall of the free surface end- 
ing in a vertical tangent at the cusp point. For clarity, 
the lighting is chosen such that the free surface appears 
opaque, so the falling jet is only indicated symbolically 
to guide the eye, but not visible directly. 

Increasing the flow rate still further, there is a second 
critical value where the stationary proflle of Fig||a ceases 
to exist and a sheet of air shoots out from the tip of the 
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cusp. The bottom picture shows this dynamical structure 
l/60th of a second after the cusp has become unstable. 
A thin air sheet now forms the wall of a transparent fluid 
cylinder. The details of this dynamical structure, such 
as the bell-shaped opening at its lower rim, is not the 
subject of this paper, but only the instability of the static 
shape leading to it. The cylinder eventually grows to 
about ten times the length shown, and is unstable to 
the formation of bubbles at its lower end, so the liquid 
pool quickly becomes contaminated by bubbles of a broad 
variety of sizes. 




FIG. 2. The local shape of the cusp, cut perpendicular to 
the sheet of air. The variable y is the distance from the tip. 

Since the air sheet near the cusp is of micron size, the 
curvature with which it is wrapped around the impinging 
jet is of no consequence and the cusp can be viewed as 
a two-dimensional object. In this spirit, Joseph et al. 

performed experiments with a two-dimensional flow 
produced by two counter-rotating cylinders submerged 
below the surface of a very viscous liquid. If the cylinders 
are placed sufRciently close to each other, a cusp forms in 
a symmetrical position between the cylinders . Letting 
the distance between the rollers going to zero, Jeong and 
Moffatt found a family of exact solutions to this problem, 
obeying the local scaling form 
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where k is the curvature at the tip (see Fig.||) . The struc- 
ture of the solution (|l|) is typical for flows near singular- 
ities |Pt|-|i3|], which involve very small scales. Physically 
(|l|) means that the shape of the interface is independent 
of scale, up to a rescaling of the axes. The other crucial 
property of singularities is that their shape is universal, 
i.e. independent of the particular type of flow that gener- 
ates the cusp. Thus our theory, based on the stability of 
such a singular structure, will be equally general. Indeed, 



Antonovskii discovered yet another class of exact so- 
lutions, but where the cusp is formed on the surface of a 
circular bubble. A local analysis reveals that the scaling 
function 

is identical to the one in except for a different value 
of the numerical constant a, confirming the expectation 
that the fiow on small scales is universal, independent of 
the particular features of the driving flow. 

In all solutions, the tip curvature k grows exponen- 
tially with the capillary number Ca — rjU where 77 
is the fluid viscosity, 7 the surface tension, and U \s a, 
typical velocity scale of the external flow. Thus as the 
strength of the driving flow increases relative to the sur- 
face tension, the size of the tip may easily reach micro- 
scopic dimensions [Q if the effect of the air is not taken 
into account. Without it, stable solutions are thus pre- 
dicted to exist for all capillary numbers, in disagreement 
with experiment. Moffatt suggested |15| that this is be- 
cause all previous analyses neglect the viscosity Ary of the 
air being drawn into the cusp by the flow Uy"-* (y) parallel 
to the cusp surface. The air entering a narrow space and 
having to escape again generates a so-called lubrication 
pressure piubiv) inside the cusp, whose derivative with 
respect to the distance y from the cusp is 



pU = 3Xvu(^\y)/h^y) 



(2) 



by Reynolds' theory pd. Since the cusp narrows like 



y 
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the lubrication pressure pushes the walls 



apart according to p;„f, ~ , just as it would keep sep- 
arated to narrowly spaced mechanical parts. 

Figure ^ proves by direct numerical simulation that 
this is enough to destroy the stationary solution found 
for A = 0. We use a boundary integral code (l^Jl^, opti- 
mized to resolve the cusp between two merging cylin- 
ders neglecting the fluid inertia. Starting from 
Antonovskii's solution with kq — 10'*, A is increased in 
steps of 2.5 • 10~^, pushing the interface forward, but only 
every fourth profile is shown. At A = 5.5 - 10^"', no more 
stationary solution is found, but instead air enters the 
fiuid forming a narrow sheet, as seen in Fig. |l] and ob- 
served in earlier experiments [ pO[ An important conse- 
quence is that in a physically correct description which 
incorporates the effect of the air (or some other gas at- 
mosphere), molecular dimensions are never reached, so 
continuum theory remains valid throughout [ pl| . 

To describe the influence of the air analytically, note 
that the extra transverse velocity fleld u"x\y) generated 
by the air pressure can simply be added to u'x\y) as 
given in [Q, since Stokes' equation is linear. Geometri- 
cally, the cusp looks like a two-dimensional crack enter- 
ing the fluid, a problem well studied in linear elasticity 
p2[ . Borrowing Muskhelishvili's result, we can now write 

v^\y) as 
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ui^Hy) = / piy')miy'/y)dy', 
Jo 

m{x) = (l/27r) ln((l + - V^)). 



(3) 




FIG. 3. A boundary integral simulation of a bubble in the 
flow proposed by Antonovskii for e = 5 and Ca = 0.0992 |^^. 
The undisturbed bubble radius is used to non-dimensionalize 
all lengths in the problem. As A is increased, the tip is pushed 
forward, but becomes narrower at a given y. The lowest pro- 
file is non-stationary. The inset shows the critical value of A 
beyond which there is no more stationary solution for a given 
curvature. 

But our free-surface problem is of course non-linear, 
since the free surface has to follow the streamlines of the 
flow, which are modified by Ux^\ Namely, the inverse 
slope of the interface is 



h' = (4°) + 



(4) 



where u^'^ and u'^^ are known 



from h as outlined above. Note that while u^'' has to 
point inward towards the cusp, wi'^'' results from the lu- 
brication pressure and points away from the cusp. Thus, 
at a given distance y from the tip, h' becomes smaller 
and the channel narrows. Owing to (|^) the lubrication 
pressure is increased, further increasing ui^^ , so this non- 
linear feedback eventually destroys the cusp solution, as 
seen in Fig. ||. 

It is extremely useful to recast equations in the 

scaling variable ^ = y^}^"^, cf. (|^). First, from (^) and 
since u^'^ is a constant up to logarithmic corrections, 
must scale like k~^/^k^/'^ = k~^/^. From (||) piub is esti- 
mated as piub ^ Ak, and thus ui"^ ^ Ak^/^ from integrat- 
ing once. The two opposing velocities become compara- 
ble at some critical value of the parameter r — Xn^^'^. 
Thus (|)-(|) can be recast in similarity variables, leading 



to an integral equation for the correction Hc{(,) to the 
unperturbed surface profile H(^): 



3rC 



dr] (5) 



where M'{rj) = m{ri). It is a simple matter to solve (|^) 
numerically, giving increasingly large corrections Hc{^) 
to the profile as r is raised. Since He is negative, the 
denominator in the integrand of (|^) decreases, leading 
to a further increase in the absolute magnitude of the 
correction, in accordance with the qualitative argument 
given above. Owing to this nonlinear feedback, a solu- 
tion ceases to exist above a critical value of r, which has a 
weak dependence on k due to the logarithmic dependence 
of u^y'^ on its argument. Hence for the flow parameters of 
Fig. ^ we predict that the cusp becomes unstable when 
the curvature reaches a critical value of Kcr ~ 0.45A^^/"^. 
This approximation is hardly distinguishable from the re- 
sult of the full solution of (^, which in the insert of Fig. 
^ is seen to be in good agreement with numerical simula- 
tions for various values of A. Because of the relationship 
between curvature and capillary number, this translates 
into the anticipated critical value Ca^r above which sta- 
tionary solutions are no longer found. At low viscosities, 
the capillary number never even reaches the critical value 
for the formation of a cusp, so an unperturbed water jet 
does not entrain air 

In conclusion, we have incorporated the effect of an 
outer fluid like air into the theoretical description of a 
cusp. This allows for the flrst quantitative description of 
air entrainmcnt through surface singularities. A descrip- 
tion of the resulting sheet of air and its stability remains 
to be done. Other future challenges include the analysis 
of its close three-dimensional relatives, namely jet forma- 
tion out of bubble tips "tip-streaming" , Taylor cones 
"electric jets" |Q, or spouts formed by planar interfaces 
"selective withdrawal" [psl . 
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